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When applied to different input states, an imperfect quantum operation yields output states with 
varying fidelities, defined as the absolute square of their overlap with the desired states. We present 
an expression for the distribution of fidelities for a class of operations applied to a general qubit 
state, and we present general expressions for the variance and input-space averaged fidelities of 
arbitrary linear maps on finite dimensional Hilbert spaces. 
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I. INTRODUCTION 

In quantum control and quantum information theory 
one attempts to control the dynamics of a quantum sys- 
tem, such that the net mapping of a quantum state by 
the dynamics yields a specific output state. If the system 
is known to initially populate a specific state, the success 
of a given operation on that state is determined by the 
square of the overlap between the final state and the de- 
sired output state, while, more generally, the successful 
implementation of a certain map, e.g., a quantum gate in 
a quantum computer, should be judged by an evaluation 
of this overlap for all possible input states. Imperfections 
may be due to a variety of reasons, such as dissipative 
coupling to auxiliary degrees of freedom and imprecise 
knowledge of the physical parameters characterizing the 
system. A wide range of methods have thus been applied 
to counter such effects: composite pulses, quantum con- 
trol, bang-bang control, error-correcting codes, and use of 
decoherence-free subspaces.T o optimize these methods 
it is necessary to have definite functionals of the gate op- 
eration, that one can determine and, hopefully, improve 
by suitable variation of accessible control parameters. 

In this paper, we address the distribution of fidelities 
obtained when the input states are taken uniformly from 
the full Hilbert space of the physical system or from a 
suitable subspace, to which physical circumstances may 
restrict the initial state. In Sec. II we present a derivation 
of the fidelity distribution of unitary gates applied to a 
qubit, i.e., a two-level quantum system, and an extension 
to the case of a unitarily diagonalizable linear map. In 
Sec. Ill, we review our recent derivation [l|] of the mean 
value of the fidelity of an arbitrary linear map, and in Sec. 
IV, we determine the variance of the fidelity distribution 
of such linear maps. Sec. V concludes the paper. 



II. FIDELITY DISTRIBUTION FOR A 
ONE-QUBIT GATE 

Consider a two-level system, a qubit, subjected to a 
unitary operator U , while the desired operation on the 



system is given by the unitary operator Uq. Under ideal 
circumstances, the unitary operator M = UqU — Uq^U 
is the identity operator, and the fidelity of the quantum 
operation \{'ip\U^U\'ip)\'^ = 1 for all normalized states ji/')- 
Under less ideal circumstances, however, M is not the 
identity, but it can be diagonalized by a unitary transfor- 
mation, and has two complex eigenvalues e"^" and e"^^ , 
with (po y^ (pi. We observe that the fidelity of the gate 
operation is unity when tp is equal to either of the eigen- 
vectors |0), |1) of A/, while expanding \tp), we get the low- 
est fidelity, / = Ke*-^" + e"^^)/2\^ = cos^ ((0i - (/)n)/2), 
for any equal weight superposition of these eigenvectors. 
In a Bloch sphere picture with the eigenvectors as the 
north and south pole, the fidelity is only a function of 
the polar angle 6 of the input state, \tp) = cos(0/2)|O) -I- 
e''^sin(6'/2)|l), f{9) = \e"''« cos^{9/2) + e"''^ sm^{d/2)\^, 
and with a uniform distribution of states over the surface 
of the Bloch sphere, and the corresponding polar angle 
distribution Pg ^ -^ sin 6 one readily determines the fi- 
delity distribution, Pf = J2i \df{0)/de\l^Pg,, where the 
contributions to the sum come from symmetric polar an- 
gles above and below the equator with the same fidelity. 
The expression for the fidelity distribution of a unitary, 
erroneous qubit gate thus reads. 
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/ - cos2 ( ^^i^ 



(1) 



for cos2 ((</.! - 0o)/2) < / < 1. 

We shall not attempt a derivation of the equivalent 
distribution of fidelities in higher Hilbert space dimen- 
sions. A related problem, dealing with the distribution 
of matrix elements of a Hermitian operator, was stud- 
ied by von Neumann 0, and illustrates how the solu- 
tion of such a problem breaks up in a large number of 
different cases. We shall, however, extend our analy- 
sis to the special situation where a quantum system is 
known to initially populate a two-dimensional subspace 
of a three-dimensional Hilbert space, including an excited 
state with suitable interaction properties, such that ex- 
citation from the lower qubit states to the excited states 



can be used to communicate between different quantum 
systems, e.g., by an excited state dipole-dipole interac- 
tion. In [3| a robust one-qubit gate scheme was proposed, 
in which both qubit states are simultaneously coupled to 
the excited state, giving rise to a dark state superposi- 
tion, |0), and a bright state superposition, |1), with de- 
structive and constructive interference of the couplings. 
Transferring the bright state into the excited state and 
back with laser fields with different phases implements a 
controllable phase in the dark/bright basis, equivalent to 
an arbitrary qubit rotation in the computational qubit 
basis [3 . If the coherent coupling of the bright and ex- 
cited state is imprecise and does not return the popula- 

/ 1 o\ 
tion fully to the low lying states, we obtain U — { o a j ] 

V 7* /3/ 

in the basis {|0), |1), |e)}. This implies that the effect on 
the qubit space of the application of U can be written in 
terms of U and the projection operator P on that space, 
PUP = (oq) in the basis {|0), |1)}. a is a complex 
number, and if \a\ — 1 the fidelity distribution is given 
by Eq. ^ with the difference between arg(Q;) and the 
desired phase shift by Uq replacing (j)i — <f)Q. 

An operation with |a| < 1 yields a special example of 
the slightly more general case, where the state vector is 
mapped, |-0) i-^ N\ip), by a normal 2x2 matrix N . The 
matrix M = UqN is also normal, i.e., it can be diagonal- 
ized by a unitary transform with eigenvalues Ao, Ai that 
we can arrange such that |Ao| < |Ai|. Employing the ex- 
pansion on the corresponding eigenstates, parametrized 
by Co = cos(0/2) and ci = e*"* sin(6l/2), where cj) S [0, 2tt] 
and 9 e [0, tt], we have / = |AoP/a, where 

/a = cos^ I + sin^ ^|A|2 + cos2 ^ sin^ ^(A -^ A*), (2) 

and A = Ai/Aq. With the same argument as above we 
need 

and we obtain 

|AopP/-||A|2-l-|l-A|2cos^f\ 
where, from p|. 



cos 6' — 



|A|2 - 1 ± 2v//A|l-A|2-Im(A)2 



II -AI- 



The fidelity is a non-monotonic function of the polar 
angle 9, and adding the contributions from different 
polar angles yielding the same fidelity, we finally obtain 
the probability distribution (recall |AoP < |Aip). We 
have to specify the results for two different cases. 



If |Ao- iAil < i|Ai|, then 
Pf 
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FIG. 1: Theoretical (line) and numerical (histogram) fidelity 
distribution for a matrix with eigenvalues Ao ~ 0.7 • e"^' * and 
Ai = 0.8 ■ e^'^'^'"'. In accordance with ((5]) the average fidelity 
is 0.28. 



and if |Ao — ^Ai| 



>5lAil 



then 



p = J |Ao-Ai| V7=7o' 



2|Ao-Ai| 77^^' 



/o < / < |AoP, 
|AoP</<|Aip, ^> 



|AoP</<|Aip, 



(3) 



where /o = Im(AoAt)V|Ao-Ai|2. 

Fig. 1. shows the fidelity distribution ^ for the case 
of Ao = 0.7 • e*^/^ Ai = 0.8 • e'^'^/^ The vertical bars 
in the figure indicate a histogram obtained by drawing 
states on the Bloch sphere at random and binning their 
individual fidelities. We observe that the results agree 
very well. 

Having the full fidelity distribution at our disposal, 
we can see how changes in the map, reflected in the pa- 
rameters {Ao,Ai}, change the fidelities, and we can op- 
timize according to different criteria: minimum value, 
mean value or some higher moment, which can readily 
be obtained from the distribution. 

In the following, we shall limit our analysis of the fi- 
delity distribution to a characterization of its mean value 
and variance. Rather than obtaining these numbers from 
integrals over the distribution ([51i)) . (/) = r^'""== Pffdf 

and (p) = r ""■"= Pfpdf, we shall provide a more gen- 
eral expression, valid for any map \'ip) i-^ N\ip), in any 
finite dimensional setting. 



III. THE AVERAGE FIDELITY MEASURE 

Derivations have been given in the literature for the 
average fidelity of qubit ; 4?| and qudit operations |l|, [^ , 
as well as for their evaluation as a sum over a properly 
chosen discrete set of states |6i, i7i] . The latter approach 
is particularly relevant in connection with quantum 
process tomography [8|, which provides a procedure for 
determining the quantum operation acting on a system 
from experimental observations. 



For any linear operator M acting on an n-dimensional 
complex Hilbert space, the uniform average of | {tp\M\tp) p 



over state vectors 
given by 



on the unit sphere S' 



2n-l 



(/> 



S2--1 



|M|V)|^dF 



-^-L_^[Tr(AfMt) + |Tr(Af)|2], 



(5) 



where dV is the normalized measure on the sphere. 

Recent detailed proofs of this result are given in [l|, M, 
llO|, and it is also readily verified using a recent result 
for the averages of general polynomials of state vector 
amplitudes over the unit sphere [U\ . 

The expression ([5]) is readily evaluated for any ma- 
trix M, and in particular, for the qubit example with 
eigenvalues Aq and Ai, we get the simple result, (/) = 
(|AoP + |Ai|2 + Re(AoA|)) /3. If M is unitary, MAft is 
the identity matrix with trace equal to the dimension n 
in Eq.(5). 

In the following we recall a few useful results, derived 
inQ. 



Subspace averaged fidelity of a unitary 
transformation 



In a number of quantum information scenarios, auxil- 
iary quantum levels are used to mediate the desired op- 
erations. In these protocols, the auxiliary levels of the 
quantum system may, with certainty, be unpopulated be- 
fore the process and, consequently, we should only aver- 
age the fidelity over the relevant input states. Since the 
final state is, ideally, also in the same subspace, we con- 
sider the matrix M = {PUIp){PUP), where P is the 
projection operator on the relevant, quantum informa- 
tion carrying subspace §, and Uq represents the desired 
unitary evolution of the relevant subspace. In a ma- 
trix notation, the outermost applications of the projec- 
tion operator P amounts to the extraction of the square 
rirei X rirei submatrix M^ei with the relevant columns and 
rows of the full matrix M , and to compute the mean fi- 
delity over the subspace, we employ (O for this reduced 
matrix: 



integrated acceptance probability over the input Hilbert 
space, ^ \\PUP\il))\\'^dV . We hence obtain the average 
conditional fidelity [20] 



Fr 



mpui 



PUP 



-m' 



\pup\'4>)\\ 



pc/p|V)ir'^" J\\PUP\ij)W^dv 

J\{il;\PUlPUP\ij)\'^dV 
J(^\PWPUP\i})dV 

1 Tt{uIpUPWPUo) + \Tt{U^PUP)\^ 



dV 



Tlrel + 1 



TiiWPUP) 



where the numerator is evaluated using ^ and the de- 
nominator follows from J {tp\M\ip)dV = Ti{M)/n. 



B. Average fidelity of a general quantum operation 

Our quantum system may not be fully isolated from its 
surroundings, and ancillary quantum systems may play 
significant roles in various quantum information process- 
ing scenarios: quantum memory protocols in a very ex- 
plicit manner involve an extra quantum system, quantum 
teleportation requires an extra entangled pair of systems, 
and in quantum computing with trapped ions, a motional 
degree of freedom is used to couple the particles. The an- 
cillary systems are ideally disentangled from the qubits 
before and after the process, but in general they act as 
an environment and cause decoherence of the quantum 
system of interest. This forces us to generalize the formal- 
ism and take into account the general theory of quantum 
operations, according to which the mean dynamics is ac- 
counted for by density matrices which transform by com- 
pletely positive maps. According to the Kraus represen- 
tation theorem, any completely positive trace-preserving 
map Q admits the representation 



gip) = J2GkpGl 



(7) 



where ^j. G\Gk — In is the n x n identity matrix [8|. 

If the pure input state p = \^p){ip\ is mapped to the 
output state g{p) the mean fidelity with which our oper- 
ation yields a unitary transformation Uq is 



F = 



-. -—[TTiMrelMli) + \TTiMrel)\^], (6) „ f / nrrt m I \ f H^TT U \ ^^r 



If population leaks to the auxiliary levels, Mrei is not 
unitary, and hence both terms of (O have nontrivial val- 
ues. 

If a measurement assures that the final state does not 
populate the complement to S, it is meaningful to de- 
fine the average conditional fidelity Fc- If we only ac- 
cept the state if it is in §, the squared overlap between 



J2 1 mMt^m'^dv 

-^{Tr(^MtA4)-fE|T>-(M.)P}, 



(8) 



the conditional, renormalized state 



PUP\-4i) 



and the 



ideal state f/o^|^) must be weighted with the accep- 
tance probability ||PC/P|V')||^ and renormalized by the 



where Mk = C/oGfe- 

Equation ([5]) enables the calculation of the average fi- 
delity of any quantum operation, as soon as it has been 
put in the Kraus form. For examples, see [l|. 



IV. THE VARIANCE OF THE FIDELITY 
DISTRIBUTION 

To further characterize the fidehty distribution, we 
shall now obtain an explicit formula for the variance of 
the fidelity distribution. Note that it is not meaningful 
to define the variance of a general quantum transforma- 
tion governed by the Kraus form, since the density matrix 
formulation already incorporates an averaging procedure, 
due to the trace over unobserved degrees of freedom of 
the surroundings of the quantum system. The density 
matrix evolution may be understood, and simulated , by 
an average over randomly evolvin g w ave functions [12l |. 
but this unravelling is not unique [l3| , and only if mea- 
surements are actually carried out on the surroundings, a 
specific stochastic dynamics of the state vectors happens, 
and it is possible to assign a definite fidelity distribution, 
and variance of this distribution, to a given process. 

We shall proceed and derive a relation for the uniform 
average of /2, i.e., |('(/'|M|V')|'^, with A/ = U^N, for linear 
state vector maps \ijj) i— s- N\ijj), where N is any linear 
operator. 

First, we show that for a Hermitian operator S, the 
uniform average of |(-!/'|S'|V^)|'* over state vectors |V') on 
the unit sphere S*^""^ in C" is given by 



I 



msm'dv 



n{n 4 

[6Tr(S'*) + 8Tr(S'3)Tr(5) 4 
-I- 6Tr(52)Tr(S')2 + Tr{S)^ 



l){n + 2){n- 
3Tr(S'2)2 



3) 



(9) 



Our demonstration of this result proceeds along the 
lines of our demonstration of Eq.(l5]), given in more detail 
in [l| . First we note the invariance of both sides of ^ 
under conjugation by any unitary operator S — *■ U^ SU, 
which allows us restrict the analysis to a diagonal ma- 
trix A with real eigenvalues Ai, . . . , A„. Letting L denote 
the left-hand side of (O, we observe that L{A) is a ho- 
mogeneous polynomial of degree 4 in the real variables 
Ai, . . . , A„, which is invariant under the exchange of any 
two Ai and Aj . It is easy to demonstrate that the set 

{Tr(A'*),Tr(A3)Tr(A),Tr(A2)2,Tr(A2)Tr(A)2,Tr(A)'^} 

spans all such polynomials i(A), and by evaluating the 
integrals 



S2n.- 



S2n.- 



S2n.- 



52,1-1 



'dV 



4|„ |4 



24 



c.l'^dV 



n{n+ l)(n + 2)(n + 3) 
4 



n{n + l){n + 2){n + 3) 
' ■'' n(n + l)(n + 2)(n + 3) 



'\c,nck\'h\'dV 



n{n + l){n + 2){n + 3) 
1 



n(n-|-l)(n-|-2)(n + 3)' 
(10) 



for different i,j,k,l, and by choosing choosing five ap- 
propriate matrices for which the integrals are readily ob- 
tained, we finally obtain the coefficients in ([9]) by solving 
a linear system of equations. 

Note that ([9]) also holds for an anti- Hermitian matrix 
A, since L{A) = L{iA) ^ R{iA) = R{A), where L and R 
denote the left- and right-hand sides of ^, respectively. 

Having obtained the uniform average of |(';/'|S'|V')|^j 
where S is Hermitian, we now proceed to the general 
case. We decompose the general matrix M = S + A 
as a sum of a Hermitian and an anti-Hermitian matrix, 
denoted by S and A, respectively, and we note that 



mM\ip)\uv = 
Msm'+mAM^+2Msm^mAM'dv. 



52„-l 



(11) 



The first two terms on the right-hand side of (|TT|) are 
readily evaluated using ^ and the explicit expressions 
S = {M + Mt)/2 and A ^ {M - M^)/2. 

To calculate the third term, we use the conjugation 
invariance to diagonalize the Hermitian part S, and thus 



msm\MAm^dv ^ 



(12) 



where A = USU^^ is a diagonal matrix with elements 
Ai, . . . , A„ e R, and A = UAU~^ is an anti-Hermitian 
matrix, which is not necessarily diagonal. Expanding the 
state vectors in the eigenbasis of S, {ip) — '^Cj\j) and 
employing the invariance of the integral (|12p under the 
transformation Cj — > exp(i0j)cj, it follows that 



MAm^mAM^dv = 

-1 

/ Y] \Caf\Ci3fc*c'^CjCkXaXpA,jAliX 

/ 02,1-1 

•^^ a,l3,i,j,k,l 

1 

(SijSki + SikSji ~ SijSki6ii)dV 



n{n + l){n + 2){n + 3)' 



{24^ A^l^p + Y^ XkXiiAuA*^ + |i. 



ijkl 



4^' [(^' + 3A^A,)|^|' + (3A2 + 2A,A,)(^i*, + |i,,f ) 



Y! [^^Aj-I^fcl' + (Afc + 4A,Afe) (A,,i*^. + |i,,f )] } 

(13) 



ijk 



where the last step follows upon insertion of the ex- 

indi- 
cates that all indices must be different. After a lengthy, 
but straightforward calculation, p3p can be rewritten in 



terms of traces of products of powers of A and A. In- 
voking trace invariance, we can replace A and A with S 
and A, respectively, and insert the explicit expressions 
for S and A in terms of M and M^ . Collecting terms 
in (|11[) . we finally obtain the man value of the squared 
fidelity, i.e., the uniform average of |(?/'|M|-0)|'* for the 
linear operator M , 



1 



where P is the corresponding projection operator on the 
relevant subspace. 



\{i;\MmUV^—— 

I n[n + l)(n + 2)[n + S) 

[ATt{M^M^^) + 2Tr{MM'<MM^) 
+ 4Tr(A//)Tr(MAf^2) ^ 4Tr(Af^)Tr(A/2Af^) 
+ Tr(A/2)Tr(Af^2) ^ 2Tr(Af A/^)^ 
+ Tr(Af2)Tr(Af^)^ + Tr{MfTT{M^^) 

+ 4Tr(A//)Tr(A/^)Tr(A/Af^) + Tr(M)2Tr(A//^)2] . 

(14) 



Combining ([5]) and (|14p . we obtain the variance of the 
fidelity distribution for a linear transformation \i/;) i~> 
N\ip), using 



^} = if) if) 

MM\i,)\UV- 



MM^'dV 



with M ~ UqN, where Uq and N are the desired and ac- 
tual evolution operators of the system. As in Sec. IIII A| 
this result is easily extended to the case where the average 
is solely performed over a subset of input states by sim- 
ply replacing n with rirei, and M with A/^ej = PUqPNP, 



SUMMARY 



We have in this paper derived an expression for the fi- 
delity distribution, applicable to a normal linear state 
vector transformation in a two dimensional complex 
Hilbert space. We have also presented simple and com- 
pact expressions for the average fidelity of a general quan- 
tum operation, and for the variance of the fidelity of a 
linear state vector map. Such simple expressions consti- 
tute a good starting point for further analysis, e.g., for 



the achievements of error correcting codes [15| , decoher- 
ence free subspaces [la, [l3l , and protection of quantum 
information by dynamical decoupling [l8|. Our expres- 
sion can also be handled and generalized analytically, as 
illustrated by our study in 1] of a ii'-qudit register, which 
provides insight into the scaling of errors. This may have 
applications in quantum error correction, the capacity 
of quantum channels, and the way that, e.g., commu- 
nication with quantum repeaters [l9| and entanglement 
distillation should optimally be carried out. Although we 
did not consider that possibility here, we note that the 
ability to restrict averages to subspaces may also enable 
generalization of our formalism to deal with non-uniform 
averages, assuming nontrivial prior probability distribu- 
tions on the Hilbert space. 
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